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ABSTRACT 



In a previous paper (Yuan 1999, hereafter Paper I), we find that the outer bound- 
2 ! ary conditions(OBCs) of an optically thin accretion flow play an important role in 

I determining the structure of the flow. Here in this paper, we further investigate the 

influence of OBC on the dynamics and radiation of the accretion flow on a more de- 
tailed level. Bremsstrahlung and synchrotron radiations amplified by Comptonization 
^ ' are taken into account and two-temperature plasma assumption is adopted. The three 

■ OBCs we adopted are the temperatures of the electrons and ions and the specific an- 

gular momentum of the accretion flow at a certain outer boundary. We investigate the 
individual role of each of the three OBCs on the dynamical structure and the emergent 
spectrum. We find that when the general parameters such as the mass accretion rate M 
and the viscous parameter a are fixed, the peak flux at various bands such as radio, IR 
and X-ray, can differ by as large as several orders of magnitude under different OBCs 
in our example. Our results indicate that OBC is both dynamically and radiatively 
important therefore should be regarded as a new "parameter" in accretion disk models. 

As an illustrative example, we further apply the above results to the compact radio 
source Sgr A* located at the center of our Galaxy. Advection-dominated accretion flow 
(ADAF) model has been turned out to be a great success to explain its luminosity 
and spectrum. However, there exists a discrepancy between the mass accretion rate 
favored by ADAF models in the literature and that favored by the three dimensional 
hydrodynamical simulation, with the former being 10-20 times smaller than the latter. 
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By seriously considering the outer boundary condition of the accretion flow, we find 

that due to the low specific angular momentum of the accretion gas, the accretion in 
Sgr A* should belong to a new accretion pattern, which is characterized by possessing a 
very large sonic radius. This accretion pattern can significantly reduce the discrepancy 
between the mass accretion rate. We argue that the accretion occurred in some detached 
binary systems, the core of nearby elliptical galaxies and active galactic nuclei(AGNs) 
very possibly belongs to this accretion pattern. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: active - 
Galaxy: center - hydrodynamics - radiation mechanisms: thermal 

1. INTRODUCTION 

The dynamics of an accretion flow around a black hole is described by a set of non-linear diff'er- 
ential equations. So this is intrinsically an initial value problem, and generally the outer boundary 
condition(OBC) has a significant infiuence on the global solution. In the standard thin disk model 
(Shakura & Sunyaev 1973), all the differential terms are neglected therefore the equations are re- 
duced to a set of algebraic equations which don't entail any boundary condition at all. However, 
in many cases it has been shown that the differential terms, such as the inertial and the horizontal 
pressure gradient terms in the radial momentum equation and the advection term in the energy 
equation, play an important role and can not be neglected (Begclman 1978; Begelman & Meier 
1982; Abramowicz et al. 1988). This is especially the case for an optically thin advection domi- 
nated accretion flow (ADAF) (Ichimaru 1977; Rees et al. 1982; Narayan &; Yi 1994; Abramowicz 
et al. 1995; Narayan, Kato &; Honma 1997; Chen, Abramowicz &; Lasota 1997). It is the existence 
of the differential terms that will make OBC an important factor to determine the behavior of an 
accretion flow. 

This expectation has been initially confirmed in our Paper I. In that paper, taking optically 
thin one- and two-temperature plasma as examples, we investigated the influence of OBC on the 
dynamics of an accretion flow. Wc adopted the temperature and the ratio of the radial velocity to 
the local sound speed (or, equivalently, the angular velocity riout) at a certain outer boundary rout 
as the outer boundary conditions and found that in both cases, the topological structure and the 
profiles of angular momentum and surface density of the flow differ greatly under different OBCs. In 
terms of the topological structure and the profile of the angular momentum, three types of solutions 
are found. When Tout is relatively low, the solution is of type I. When Tout is relatively high and 
the angular velocity ilout is high, the solution is of type II. Both types I and II possess small sonic 
radii, but their topological structures and angular momentum profiles are different. When Tout is 
relatively high and the angular velocity is lower than a critical value, the solution is of type III, 
characterized by possessing a much larger sonic radius. For a one-temperature plasma or ions in a 
two-temperature plasma, the discrepancy of the temperature in the solutions with different OBCs 
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lessen rapidly away from the outer boundary, but the discrepancy of the electron temperature at 
''out persists throughout the disk. This is because when the accretion rate is low, the electrons 
are basically adiabatic, i.e., both the local radiative cooling and the energy transfer from the ions 
to the electrons can be neglected compared with the energy advection of the electrons. While 
for a one-temperature plasma or ions, the local viscous dissipation in the energy balance plays an 
important role, thus their temperature is more locally determined. 

In the present paper we will focus on the two-temperature accretion flows. Our first aim is to 
extend our study in Paper I by including synchrotron radiation as well. This is a very strong local 
radiative cooling mechanism in the inner region of a disk when a magnetic field is present. We will 
check whether our conclusion of Paper I still holds or not in this case. 

The main aim of this paper is to calculate the emergent spectrum under different OBCs and 
to investigate the applications. Wc want to probe whether and how the spectrum is dependent 
upon the OBCs (Sections 2 & 3). We find that each of the three OBCs has a significant influence 
on the emergent spectrum. This result indicates the importance of considering the initial physical 
state of the accretion flow. In Section 4 we apply our theory of OBC to Sgr A*. We find that 
the discrepancy between the mass accretion rate favored by the numerical simulation and that 
required in all ADAF models in the literature is naturally solved by seriously considering the OBCs 
of accretion flows. The last section is devoted to summary and discussion, where the promising 
applications in nearby elliptical galaxies and Active Galactic Nuclei (AGNs) are shortly discussed. 



2. MODEL 

Wc consider a steady axisymmetric accretion flow around a Schwarzschild black hole of mass 
M. Paczyhski &: Wiita (1980) potential = —GM/{r — Vg) is used to mimic the geometry of the 
hole, where Vg = 2GM/(? is the Schwarzschild radius. The standard cc-viscosity prescription is 
adopted. We assume that all the viscous dissipated energy is transfered to ions and the energy 
transfer from ions to electrons is provided solely by Coulomb collisions. So the plasma has a two- 
temperature structure in the present optically thin case. A randomly oriented magnetic field is 
assumed to exist in the accretion flow. The total pressure is then taken to be 

P = Pgas +Pmag, (1) 

where pgas is the gas pressure and Pmag is the magnetic pressure. For simplicity, we assume the ratio 
of gas pressure to the total pressure /3 to be a global parameter independent of radius r. Under the 
optically thin assumption, the equation of state can be written as: 

Pgas = Pp=Pi+Pe = — —Ti + — —Te, (2) 

Hi rUfj, He rUfj, 

and 

Pmag = (1-/3)P= g^- (3) 
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Here and hereafter subscripts i and e indicate the quantities for ions and electrons, respectively. 
The mean molecular weight for ions and electrons are: 

^li = 1.23 = 1.13, (4) 

respectively. We take the shear stress to be simply proportional to the pressure, i.e., shear stress=- 
ap. The hydrostatic balance in the vertical direction is also assumed. 

Under the above assumptions, the set of height-integrated equations describing the behavior 
of accretion flows read as follows. 

- AnrHpv = M, with H = Cg/^k = \/p/p/^k, (5) 

v^^ = -nir + n\-'yP (6) 

dr p dr 

t;(J7r^ — j) = ar-, (7) 

fdsi d fl\\ , dn 

P"^ {-r + Pi~r [~ I =Q -Qie = -apr- qie, (8) 

\dr dr \p / / dr 

P^i'^+Pe^i-X] =(lie-q~- (9) 



dr dr \p, 

All above quantities have their popular meanings. The difference between eq. (7) above and eq. 
(2.11) in Narayan, Kato, & Honma (1997) or eq. (2.13) in Narayan, Mahadevan, & Quataert (1998) 
is due to the difference of the adopted viscosity prescription (see, e.g., Abramowicz et al. 1988; 
Nakamura et al. 1997; Manmoto, Mineshige, & Kusunose 1997). We adopt this kind of viscosity 
prescription because in this case the no-torque condition at the hole horizon is automatically satis- 
fied therefore the calculation can be simplified (Abramowicz et al. 1988; Narayan, Kato, &: Honma 
1997). qie denotes the energy transfer rate from ions to electrons by Coulomb collisions which takes 
the form (Dermer, Liang & Canfield 1991): 

3 me iif'^ + {Ge + eif''^ 
lie = -— nenjO-TclnA {kTi - kTe) ^ -rp^ , (10) 

where InA = 20 is the Coulomb logarithm, and 9i = kTi/rriiC^ and 9e = kTe/rrieC^ is the dimension- 
less ions and electrons temperatures, respectively. This formula is accurate enough for our purpose 
because Qie is only a small fraction in the energy balance of both ions (eq. 8) and electrons (eq. 9) 
(Nakamura et al. 1997). Si and Se denote the internal energies for ions and electrons per unit mass 
of the gas, respectively. Since in ADAFs the ions never become relativistic while the electrons are 
transrelativistic, following Esin et al. (1997b), we adopt the following formsQ: 

1 kiTii 3 JxiTii 



7i - 1 piirin 2 piTUf, 



(11) 



^ based upon the argument of Esin et al. (1997a), the adiabatic index "fue) in Esin et al.(1997b) include the 
contribution of the magnetic density as well. But Quataert & Narayan (1999a) found that this is incorrect if MHD 
adequately describes the accretion flow. So we exclude it in the present paper. See Quataert & Narayan (1999b). 
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7e - 1 fJ-emf, Hem^ 
The expression for coefficient a(Te) is (Chandrasckhar 1939): 



a(re) -. (12) 



a(0e) = ^ 



3i^3(l/ee)+i^l(l/^_^' 



(13) 



4i^2(l/^e) 

where ^1,7^2, and are modified Bessel functions of the second kind of order 1, 2, and 3, 
respectively. 

As for the calculation of the radiative cooling q~ ^ we follow the procedure in Manmoto, Mi- 
neshige h Kusunose (1997). The considered radiative mechanisms include bremsstrahlung, syn- 
chrotron radiation and Comptonization of soft photons. Assuming the disk is isothermal in the 
vertical direction, the spectrum of unscattered photons at a given radius is calculated by solving 
the radiative transfer equation in the vertical direction of the disk basing upon the two-stream 
approximation (Rybicki h, Lightman 1979). The result is (Manmoto, Mineshige &: Kusunose 1997): 

= -^S,[l - exp(-2V3<], (14) 

where r* = {ir^/^ /2)ki,H is the optical depth for absorption of the accretion flow in the vertical 
direction with Ki, = Xu/i'^'^Bi,) being the absorption coefficient, where Xv = X;/,brems + Xi^,synch 
is the emissivity, and X;/,brcms ^ind Xj^,synch are the bremsstrahlung and synchrotron emissivities, 
respectively. Then the local radiative cooling rate q~ reads as follows: 

1~ = ^j dM'^)^F,, (15) 

where 77 is the energy enhancement factor first introduced by Dermer, Liang &; Canfield (1991) and 
modified by Esin et al. (1996) (See Manmoto, Mineshige &; Kusunose 1997 for the exact formula). 

Given the values of parameters M, M, a and /?, we numerically solve the above set of equations 
describing the radiation hydrodjmamics of a two-temperature accretion flow around a black hole. 
The equations are reduced to a set of differential equations with three variables v, Ti and Tg. The 
global solutions must satisfy simultaneously the no-torque condition at the horizon, the sonic point 
condition at a sonic radius rg and three outer boundary conditions given at a certain outer boundary 
Tout- The numerical method we adopted is the same as in Nakamura et al.(1997). 

We adopt the same procedure as in Manmoto et al. (1997) to calculate the emergent spec- 
trum. The spectrum of unscattered photons is calculated from equation (14) and the Compton 
scattered spectrum is calculated by using the formula given by Coppi k, Blandford (1990). Since 
the Comptonization is mainly occurred in the inner region of the disk r ^ lOrg, where H/r ^ 0.4, 
our "local" instead of "global" treatment of Comptonization will not cause a serious error. At 
last, the redshift due to gravity and gas motion are also considered in the following simple way 
(Manmoto, Mineshige &: Kusunose 1997). The effect of gravitational redshift is included by taking 
the ratio of the the observed energy of a photon to its energy emitted at radius r to be (1 — r^/r)^/^. 
As for the redshift due to the relativistic radial motion, we concentrate on the face-on case and 
take the rate of energy change to be 1/[1 — (u/c)^]. 
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3. RESULTS 

Throughout this paper, we set M = IO^Mq and m = M/MEdd = 10"^, where Msdd = 
lOLEdd/c^ = 2.2 X lO^^Myr"-*^ is defined as the Eddington accretion rate. Other parameters are 
a = 0.1,/? = 0.9 and rout = lO^rg. At rout) the three outer boundary conditions we imposed 
are Tj = Tout.ijTe = Tout.e and A(= v/cg ^ vj = Aout(see Paper I). This set of boundary 

conditions are equivalent to (Tout,*, Tout, f^out) according to eq. (7) because j is the eigenvalue of 
the problem rather than a free parameter. We assign (Tout,i, T'out,e5 Aout) to different sets of values 
and investigate their effects on the structure of the accretion flow and the emergent spectrum. The 
results are as follows. 

We find that the results are qualitatively the same as in Paper I although the synchrotron 
radiation, a strong local cooling term, is included in the electron energy equation. This is because 
the differential terms, which stand for the global character of the equations, still play a dominated 
role. We find that only when OBCs are within a certain range do the global solutions exist. In 
terms of the ion temperature, the range is Tout,i ~ (0.01 — l)Tvii.iai ( here T^iriai = 2GMm^/3kr 
denotes the virial temperature). The electrons temperature is slightly lower than Tout.j- As for 
the range of Oout; it must satisfy Oout ^ 0.8r2K('^out)) otherwise the viscous heating term takes 
a negative value under our viscous description. This result was first pointed out by Manmoto, 
Mineshige & Kusunose (1997). The range of OBCs slightly varies under different parameters rh,a 
and P, and it is also the function of the value of rout; as we describe in the subsequent section of 
this paper. The structures of the solutions with different OBCs are greatly different. Three types 
of solutions arc also found in terms of their topological structures and angular momentum profiles 
(cf. Paper I). Generally, when Tout is relatively low, the solution is of type I. When Tout is high 
and Aout is small (or equivalently Oout is large), the solution is of type II. Types I and II both have 
small values of sonic radii rg ^ lOrg, but their angular momentum profiles are different (see Figure 
1 of Paper I and Figure 1 in the present paper). When Tout is high and Aout is large (or equivalently 
r^out is small), the solution becomes type III, which has much larger rg. 

Figure 1 shows the effect of different Tout.i on the solution structures. The values of rout,e 
and Aout are the same, with rout,e = 1-2 x 10^ and Aout = 0.2. The solid line (belonging to 
type I solution) is for Tout,i = 2 x lO^K, the dotted line (type I) for Tout,i = 6 x lO^K, the 
short-dashed line (type II) for Tout,i = 2 x lO^K and the long-dashed line (type III) for Tout,i = 
3.2 X lO^K. Six plots in the figure show the radial variations of the Mach number (defined as 
'"/Vii^li ~ 1) + 2(7i — l)Q;^)/(7i + l)/cs so that the Mach number equals 1 at the sonic point), 
the electrons and ions temperatures, profiles of the specific angular momentum and surface density 
S (= 2pH), and the ratio of the vertical scale height of the disk H(r) to the radius, respectively. 
From the figure, we find that the discrepancy in the ion temperature is rapidly reduced as the 
radius decreases and Tj converges to the virial temperature. This is because the local viscous 
dissipation in the ions energy equation plays an important role therefore the temperature of ions 
is mainly "locally" determined, like the case of a one-temperature plasma (Paper I). However, the 
discrepancy in the surface density of the disk persists throughout the disk. 
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Such discrepancy in the solution structures results in the discrepancy in the emission spectrum, 
as Figure 2 shows. The discrepancy in the luminosity at certain individual wave band can be 
understood by referring to Figure 1. In the radio-submillimeter band, the power is principally 
due to synchrotron emission for which we approximately have Ly oc p^l^T^^'^ for the "general" 
frequency and V'pL^,^ oc (T'^'^TI for the peak frequency (Mahadevan 1997). Synchrotron radiation 
mainly comes from the inner region of the disk, <^ lOr^. So the long-dashed line in the figure 
possesses the lowest radio power is mainly due to its lowest surface density, since is almost 
the same for different solutions. For such a low mass accretion rate system as considered in this 
paper, the low frequency part of the submillimctcr to hard X-ray spectrum is mainly contributed by 
the Comptonization of synchrotron photons, and it is bremsstrahlung emission that is responsible 
for the high frequency part. The long-dashed line still possesses the lowest power in the low 
frequency band because its corresponding amount of synchrotron soft photons is the least and its 
corresponding Compton y-paramctcr y = TeskT^./rrieC^ is the smallest. In the X-ray band, the solid 
line has the highest power. This is because for bremsstrahlung emission we approximately have 
-^^brcms Tg ^/^e~'*'^/'^^'=/j^, i.e., the density instead of Tg is the dominated factor, and the solid line 
possesses the highest density among the four solutions. 

The effect of Tout,e on the structures of the global solutions is shown in Figure 3. All the three 
lines in the figure have Tout,i = 2 x IQ^K and Aout = 0.2. The solid line is for Tout,e = 1-2 x W^K, 
the dotted line for Tout.e = 8 x lO^K and the dashed line for Tout.e = 1-2 x 10^. Different with Tj, the 
discrepancy in Tg persists throughout the disk rather than converges as the radius decreases. This 
is because the electrons are essentially adiabatic for the present low accretion rate case, i.e., both 
and q~ are very small compared with the advcction term on the left-hand side of equation (8), so the 
electron temperature is globally determined. While for the ion, the local viscous dissipation in the 
energy equation plays an important role; thus, the temperature of ion is determined more locally 
than the electron. Such discrepancy in the electron temperature produces significant discrepancy 
in the emission spectrum especially in the radio, submillimeter and IR bands, as Figure 4 shows. 

The effect of modifying Aout (i-e., modifying r^out) on the structure of the accretion flows is 
shown in Figure 5. Each of the three lines in the figure corresponds to Tout,i = 2 x lO^K and 
To„t,e = 1-2 X lO'^K. The values of Aout are 0.18 (solid line, type II), 0.22 (dotted line, type II) 
and 0.26 (dashed line, type III), and their corresponding angular velocities are 0.45$7k) 0.37ilk and 
0.29rik) respectively. 

When fiout decreases across a certain critical value, ~ 0.36rik(''out) in the present case, an 
accretion pattern characterized by possessing a much larger sonic radius appears (referred to as 
type III solution). As we argue later in this paper, this accretion pattern is of particular interest 
to us. In all previous works on viscous accretion flows, only the solutions with small sonic radii ~ 
several rg have been found (Muchotreb &; Paczyhski 1982; Matsumoto et al. 1984; Abramowicz et 
al. 1988; Chen &; Taam 1993; Narayan, Kato &; Honma 1997; Chen, Abramowicz & Lasota 1997; 
Nakamura et al. 1997); therefore, this pattern is new. The reason why previous authors did not 
find this pattern is that they generally set the standard thin disk solutions as the OBCs (Muchotreb 
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&; Paczynski 1982; Matsumoto et al. 1984; Abramowicz et al. 1988; Chen k. Taam 1993; Narayan, 
Kato & Honma 1997; Chen, Abramowicz & Lasota 1997), where the specific angular momentum 
of the flow is Keplerian; or, a specific angular momentum with a value near the Keplerian one is 
set as the OBC (Nakamura et al. 1997; Manmoto, Mineshige & Kusunose 1997). 

Such transition between the sonic radii happened when f^out passes across a critical value is 
clearly shown in Figure 6. The value of the critical ilout varies with the parameters such as a and 
and it is also the functions of Tout,i, Tovx,e and rout- Generally, it decreases with the decreasing a 
and/or the increasing rout- Similar transition has been found previously in the context of adiabatic 
(inviscid) accretion flow by Abramowicz k, Zurek (1981). In that case, they found that when the 
specific angular momentum of the flow I, a constant of motion, decreased across a critical value 
lc{E) (here the specific energy E of the flow is another constant of motion), a transition from a 
disk-like accretion pattern to a Bondi-like one would happen (Abramowicz k. Zurek 1981; Lu & 
Abramowicz 1988). Here in this paper we find that this transition still exist when the flow becomes 
viscous, confirming the prediction of Abramowicz & Zurek (1981). 

The physical reason of the transition between the two patterns is obvious. When the specific 
angular momentum of the gas is low, the centrifugal force can be neglected and the gravitational 
force plays a dominated role in the radial momentum equation (eq. 6), therefore the gas becomes 
supersonic far before the horizon like in Bondi accretion. When the specific angular momentum 
is high, however, since the centrifugal force becomes stronger, the gas becomes transonic only 
after passing through a sonic point near the horizon with the help of the general relativistic effect 
(Abramowicz & Zurek 1981; Shapiro & Teukolsky 1984). 

Figure 7 shows the corresponding spectra of the solutions presented in Figure 5. The accretion 
flow belonging to the new accretion pattern, which is denoted by the dashed line, emits the lowest 
X-ray luminosity. This is because this accretion pattern possesses the largest sonic radius therefore 
the corresponding density of the accretion flow is the lowest. 

All above results are obtained for a fixed outer boundary rout = lO^rg. We also investigate 
the infiuence of increasing the value of rout on the global solutions. We find that the general 
feature are quantitatively the same. One remarkable difference is that the ranges of Tout.i and 
^out.e within which we can obtain a global solution lessen with the increasing radii. For example, 
if we increase rout from lO'^rg to lO'^rg, the temperature range will lessen from ~ (0.01 — l)Tvirial to 
^ (0.1 — l)Tviriai- If Tout is taken to be large enough, Tout,i(~ T^utje) almost be unique. However, 
the feasible range of Oout is almost constant, no matter how large rout is- Moreover, when the 
value of rout becomes larger, the transition between the accretion patterns, happened when f^out 
pass across the critical value, becomes more "obvious" in the sense that the small sonic radius 
becomes smaller and the large sonic radius becomes larger. Consequently, the discrepancy between 
the surface density of the accretion disk of the two accretion patterns becomes larger compared 
with the case of small rout, and this will further result in the increase of the discrepancy of the 
X-ray luminosity emitted by the accretion flow. This is the crucial factor to solve the puzzle of the 
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mass accretion rate of Sgr A* . 

4. APPLICATION TO Sgr A* 

Kinematic measurements suggest that the energetic radio source Sgr A* located at the center 
our Galaxy is a supermassive compact object with a mass ~ 2 — 3 x 10^ Mq. This is widely believed 
to be a black hole (Mezger, Duschl & Zylka 1996). On the other hand, observations of gas outflows 
near Sgr A* indicate the existence of a hypersonic stellar wind coming from the cluster of stars 
within several arc-seconds from Sgr A*. The wind should be accreting onto the black hole (Melia 
1992). If the flow past Sgr A* is uniform, then the mass accretion rate can be simply obtained by 
the classical Bondi-Hoyle scenario (Bondi & Hoyle 1944) as follows. Since the wind is significantly 
hypersonic, as such, a standing bow-shock is inevitable. This shock is located roughly where flow 
elements' potential energy equals kinetic energy (Shapiro 1973; Melia 1992), R\ = 2GMbh/^^) 
where Mbh is the mass of the black hole and is the wind velocity. From here the shocked 
gas is assumed to accret into the hole. Since the mass flux in the gas at large radii is mpnyjVyj, 
where nip is proton mass and Uw is the number density of the gas, the mass accretion rate is 
Mbh = -KR^mpriwVyj. 

This result is only valid for an uniform source. In reality the flow past Sgr A* comes from 
multiple sources. In this case, the wind-wind shocks dissipate some of the bulk kinetic energy and 
lead to a higher capture rate for the gas (Coker & Melia 1997). The exact value of accretion rate M 
depends on the stellar spatial distribution and can only be obtained by numerical simulation. The 
three-dimensional hydrodynamical simulation of Coker & Melia (1997) gives the average value M = 
1.6Mbh for two extreme spatial distributions, spherical and planar ones. With the available data at 
that time, the accretion rate obtained in Coker k. Melia (1997) is M ~ 10^^ gs"-*^. Considering that 
recent work suggests the Galactic Centre wind is dominated by a few hot stars with higher wind 
velocities of ^ 1000km.s^-'^(Najarro et al. 1997) rather than 700 km s^^ taken by Coker & Melia 
(1997), the more reliable accretion rate should be M (700/1000)^ x lO^^gs"! ~ 9 x 10-^MEdd- 

Assuming the accretion is via a standard thin disk, the accretion rate required to model the 
luminosity is more than three orders of magnitude lower than this value. ADAF model has been 
turned out to be a signiflcant success to model its low luminosity (Narayan, Yi &: Mahadevan 1995; 
Manmoto, Mineshige &; Kusunose 1997; Narayan et al. 1998). However, there exists a discrepancy 
between the accretion rate favored by all ADAF models in the literature and that favored by the 
simulation mentioned above, with the former being 10-20 times smaller than that favored by the 
latter (Coker & Melia 1997; Quataert & Narayan 1999b). The most up-to-date calculation gives 
M = 6.8 X lO^^MEdd (Quataert & Narayan 1999b) and this value approaches the lower limit 
considered plausible from Bondi capture (Quataert, Narayan, k. Reid 1999). In an ADAF model, 
the mass accretion rate is determined by fitting the theoretical X-ray flux to the observation. If a 
larger accretion rate were adopted, then bremsstrahlung radiation would yield an X-ray flux well 
above the observational limits. Even though we assume that signiflcant accretion mass may be 
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lost to a wind, detailed calculation shows that since the bremsstrahlung radiation comes from large 
radii in the accretion flow, the discrepancy can not be alleviated no matter how strong the winds 
are (Quataert & Narayan 1999b). 

The present research on the role of OBC tells us that we should seriously consider the physical 
state of accretion flows at the outer boundary. However, we note that in all present ADAF models 
of Sgr A*, the outer boundary condition is roughly treated. In our view, this might be an origin of 
the discrepancy between the mass accretion rate. 

Basing upon the above consideration, wc recalculate the spectrum of Sgr A*. We set the outer 
boundary at the "accretion radius" Ra = 2GM/v'^ ~ 1.5 x lO^rg. In the present case of Sgr 
A*, the temperatures of ions and electrons in the flow just after the bow-shock should equal the 
virial temperature ~ 10^^/(r/rg)if. But the value of the speciflc angular momentum at the outer 
boundary is not certain. We set it as various values and find that, when it pass across a critical 
value O.lGr^K! the transition of the accretion pattern occurs. Figure 8 shows the Mach number 
and the surface density of the two accretion patterns. The parameters of both the solid and the 
dashed lines in the figure are: Mbh = 2.5 x IO^Mq, M = 4 x IQ-'^MEdd, a = 0.1 and jS = 0.9. 
At the outer boundary Ra, the two lines possess identical temperature of 8 x lO^K, but their 
specific angular momenta are different. The solid line, which stands for our new accretion pattern, 
corresponds to relatively low angular velocity, ~ 0.15Qk, while the angular velocity possessed by 
the dashed line, which we draw for comparison, is much higher, ~ 0.46r2K- Their sonic radii are 
~ 4rg (dashed line), and ~ 6000rg (solid line), respectively. The difference of the sonic radius 
results in the discrepancy in the surface density, which further results in the difference of the X-ray 
luminosity, as shown by Figure 9. 

Figure 9 shows our calculated X-ray spectra together with the observation of Sgr A*. Here we 
assume that the bremsstrahlung radiation is the only contributor to this waveband. This assump- 
tion requires that the synchrotron radiation can not be too strong, otherwise the contribution from 
the Comptonization of the synchrotron photons will exceed that from bremsstrahlung radiation. 
This requirement can be satisfied according to the radio observation of Sgr A*. Wc will discuss this 
question in the following part. From the figure, we find that accretion rate as high as 4 x 10~*^MEdd 
is acceptable, if only the angular momentum at the outer boundary is relatively low. Compared 
with the value of 6.8 x lO^^Msdd in Quataert &; Narayan (1999b), this value is much closer to that 
favored by the numerical simulation. If the angular momentum of the flow at Ra is relatively high, 
however, as shown by the dashed line, the X-ray flux is well above the observation. 

Now the crucial moment is whether the value of the angular velocity at Ra is really low. In 

this context, we note that the three-dimensional numerical simulation by Cokcr & Mclia (1997) 
indicates that the accreted angular velocity in Sgr A* is very low, ^lout ~ 0.ir2K(?'out) for their run 
1 and « 0.2OK(?'out) for their run 2. The consistency with our result is satisfactory. 
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5. SUMMARY AND DISCUSSION 

In this paper, we numerically solve the radiation hydrodynamic equations describing an op- 
tically thin advection dominated accretion flow and calculate its emergent spectrum. We fix all 
the parameters such as a,P,M and M, and set the outer boundary condition as Tout,?, ^out.e and 
Aout(= v/Cf^) at certain outer boundary rout- Our primary focus is to investigate the effects of the 
outer boundary conditions on the structure of the global solution and, especially, on the emergent 
spectrum. 

We find that OBCs must lie in a certain range, otherwise we can't find the corresponding 
global solutions. However, this range is large enough in the sense that both the structure of the 
solution and its corresponding spectrum differ greatly from each other under different OBCs. Three 
types of solutions are also found under various OBCs, qualitatively agree with the result of Paper 
I, where magnetic field is neglected and the description to radiation is crude. The value of the sonic 
radius in type III solution is much larger than types I and II. We examine the individual infiuence 
of Tout, ij Tout, e and Aout on the structures and spectrum and find that each factor especially the 
temperature, plays a significant role. The peak fiux in the radio, IR and X-ray bands can differ 
by nearly two orders of magnitude or more in our example, although all the other parameters are 
exactly the same. 

We also investigate the effect of modifying the value of rout • We find that the feasible ranges of 
Tout,! and Tout,e lessen with the increasing radii. If the value of rout is too large, the values of Tout,i 
and rout,e are almost unique. This result seems to mean the reduction of the effect of OBC on the 
solution. However, on the one hand we expect that in some cases, such as some binary system or 
the system where a thin disk- AD AF transition is expected to occur (e.g., X-ray binaries A0620-00 
and V404 Cyg, and low luminosity AGN NGC 4258; see Narayan, Mahadevan, & Quatacrt 1998 
and references therein), the value of rout might be small. On the other hand, the feasible range of 
r^out is almost constant no matter how large rout is- 

The reason why previous works on ADAF global solution (e.g., Narayan, Kato, &: Honma 
1997; Chen, Abramowicz, Sz Lasota 1997) did not find the obvious effect of OBC on the global 
solution is that, both Narayan, Kato &; Honma (1997) and Chen, Abramowicz & Lasota (1997) 
concentrated on the dynamics of a one-temperature plasma, where the local viscous dissipation 
in the energy equation plays an important role. Due to this reason, their global solutions are in 
principle "locally" rather than "globally" determined, the effect of OBC weakens rapidly with the 
decreasing radii and the solutions converged rapidly over a small radial extent away from the outer 
boundary. We also obtained similar results with theirs for the one-temperature plasma in our Paper 
I. However, according to our result, the variation of ftout across a certain critical value will produce 
obvious OBC-dependent behavior such as the transition of the sonic radius, no matter what type 
the plasma is, one- or two-temperature. They did not find this result might be due to the fact that 
r^out adopted by them, f2K(''out) for Kcplcrian disk outer boundary condition or 0.34ilK(?'out) for 
self-similar solution outer boundary condition, is always larger than the critical value, so the effect 
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of OBC is very small hence is hard to find. 

The present study concentrates on the low-M case where the differential terms in the equation 
such as the energy advection play an important role therefore the effect of OBC are most obvious. 
When M becomes higher, the role of the local radiation loss terms in the energy balance will 
become more important. In this case, we expect that the discrepancy due to OBC in the profile of 
the temperature will lessen. However, from our calculation to the one-temperature accretion flow 
whose temperature is also principally determined locally (Paper I), and Figure 5 in the present 
paper, we expect that the flow should still present OBC-dependent behavior in, e.g., the angular 
momentum and the Mach number profiles. 

We do not include winds in the present study. Since the Bernoulli parameter of ADAF is 
positive therefore the gas can in principle escape to infinity with positive energy (Narayan &; Yi 
1994, 1995). Blandford & Begelman (1999) recently suggested that mass loss through winds might 
be dynamically important. The effect of winds on the spectrum of ADAF has been investigated 
by Quataert & Narayan (1999b). It is interesting to investigate the effect of OBC on ADAF with 
winds. We expect the result is probably similar with the present results since in that case the 
differential terms in the equation still play an important role. 

It is a meaningful problem whether a standing shock occurs in an accretion flow. Although 
some authors have set up the shock-included global solutions (e.g. Chakrabarti 1996), the result 
is generally thought not to be so convincing because in their numerical procedure Vs and j are 
treated as two free parameters instead of the eigenvalues of the problem. A necessary condition 
for the shock formation is the existence of the global solution with a large sonic radius outside the 
centrifugal barrier. This is the key of the problem. According to our result, this kind of large- 
sonic-radius solution, belonging to our type III, can only be realized when the specific angular 
momentum of the accretion flow is low. We note that this requirement to the angular momentum 
is exactly what Narayan, Kato &; Honma(1997) anticipated, and it has recently been confirmed by 
the numerical simulation by Igumenshchev, lUarionov Sz Abramowicz (1999). 

The large-sonic radius solution (type III) is a new accretion pattern since in all previous studies 
on viscous accretion onto black holes the sonic radii arc small. We find that the discrepancy in 
the mass accretion rate of Sgr A* between the value favored by the previous ADAF models in the 
literature and that favored by the hydrodynamical numerical simulation is significantly reduced if 
the accretion is via this pattern. This result hints us that such low angular momentum accretion 
may be very common in the universe. One example is the detached binary system, where the 
accretion material is the stellar wind from the companion therefore the angular momentum of the 
accreting gas is very low (Illarionov & Sunyaev 1975). Another example is the cores of nearby giant 
elliptical galaxies, where the angular momentum of the hot (accretion) gas is again assumed to be 
very small. 

We are grateful to the referee for many helpful comments and suggestions which enabled us 
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Fig. 1. — The structures of the accretion flow with different Tout,i. The sohd hne (type 1 solution) 
is for Tout,i = 2 X lO^K, the dotted Hne (type I) for Tout,i = 6 x lO^K, the dashed line (type II) 
for Tout,i = 2 X lO^K and the long-dashed line (type III) for Tout.i = 3.2 x lO^K. Other OBCs are 
^out.e = 1-2 X W^K and Aout = 0.2. The outer boundary is set at rout = lO^rg. Other parameters 
are a = 0.1, P = 0.9, M = lO^M© and M = IQ-^MEdd- The units of S and T are g cm-^ and K 
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Fig. 2.— 



The corresponding spectra of the solutions shown in Figure 1. 
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Fig. 3. — The structures of the accretion flows with different Tout.e- The sohd, dotted and the 
dashed hnes are for Tout e = 1-2 x 10^ K, 8 x 10^ K and 1.2 x 10^ K, respectively. Other OBCs are 



Touts = 2 X lO^K and Aout = 0.2. The outer boundary is set at rout 
and the units are the same as those in Figure 1. 
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Fig. 4.— 



The corresponding spectra of the solutions shown in Figure 3. 







9 







1 



2 



3 







1 



2 



3 



log(R/Rg) 



log(R/Rg) 



Fig. 5. — The structures of the accretion flows with different Aout- The sohd, dotted and the 
dashed hnes are for Aout = 0.18,0.22 and 0.26, respectively. The corresponding angular velocities 
are 0.447Qk, 0.37Ok and 0.289J7k. Other OBCs are Tout.i = 2 x IQ^K and Tout,e = 1-2 x lO^i^. 
The outer boundary is set at rout = lO^rg. Other parameters and the units are the same as those 
in Figure 1. 
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Fig. 6. — The variation of the value of the sonic radii with the angular velocity at the outer 
boundary. A transition is clearly shown. 
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Fig. 7.— 



The corresponding spectra of the solutions shown in Figure 5. 
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Fig. 8. — The radial variations of the surface density and the Mach number of the accretion flows 
with different specific angular momenta at the outer boundary i?A- The solid line, which stands 
for our new accretion pattern, is for Oout = 0.15Ok while the dashed line is for Qout = OAGQ^mK- 
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Fig. 9. — The X-ray spectrum of Sgr A*. The measured fluxes were converted to luminosity 
assuming a distance of 8.5 kpc to the Galactic Centre. The observational data are compiled by 
Narayan et al. (1998). The spectra represented by the solid and the dashed lines are produced by 
the accretion flows denoted by the same style of lines as in Figure 8. Due to the difference of the 
angular momentum of the flow at the outer boundary, the X-ray flux differs by a factor ~ 8. 



